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Abstract 

We consider three dimensional quantum electrodynamics (cQED 3 ) with massless relativistic 
fermions coupled to a compact gauge held using a combined perturbative variational approach. 
Coupling to matter renders the bare interaction between magnetic monopoles logarithmic at large 
distances, suggesting the possibility of a confinement-deconfinement transition of the Berezinskii- 
Kosterlitz-Thouless type in the theory. Our self-consistent calculation suggests, however, that 
screening effects always destabilise the confined phase, in agreement with the previous renormali- 
sation group study of the same model. 
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I. INTRODUCTION 



Compact U(l) gauge theories in three (d = 3) dimensions have long been of interest 
in high energy and condensed matter physics. In particle physics they serve as relatively 
simple models exhibiting non-perturbative phenomena such as chiral symmetry breaking and 
confinement [lL 13, [^] , believed to be crucial to our understanding of more realistic theories 
like quantum chromodynamics. In condensed matter physics the theories with compact U(l) 
gauge fields coupled to matter arise frequently in descriptions of strongly correlated electron 
systems |4] . In this case the three dimensional models are of direct significance to condensed 
matter systems in two (d = 2) spatial dimensions and at zero temperature (T = 0). 

A crucial issue in all compact U(l) theories is the confinement of 'charge' due to the 
unbinding of magnetic monopoles, which are invariably introduced by the compact nature 
of the gauge field. In a pioneering work, Polyakov |l( showed that in pure compact quantum 
electrodynamics without matter in d = 3 confinement is permanent for all values of the gauge 
coupling. The situation where the gauge field is coupled to matter is more subtle, and a 
subject of current debate. It has been argued that coupling to relativistic massless fermions 
transforms the usual Coulombic interaction between monopoles into the much longer-ranged 
logarithmic interaction at large distances □ DOB . When applied to a single monopole- 
antimonopole pair, this would suggest that monopoles may bind into dipoles, in analogy 
with the celebrated Berezinskii jjj], Kosterlitz and Thouless nj] (BKT) transition in two 
dimensions. However, while the effects of a finite density of monopoles on the BKT transition 



in d = 2 are well understood 
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lri| . the situation in d = 3 appears less clear 3|. The 
difficulty lies in the fact that while the screening in the dipole phase in d = 2 just amounts to 
renormalisation of the dielectric constant, in d = 3 it changes the form of the interaction ll3L 
UHQ. In a recent paper, two of us Q — an electrostatic ardent and a 
renormalisation group calculation to show that the interaction between distant monopoles in 
the presence of other dipoles is screened back into the Coulomb potential. Together with the 
generalisation to the case of coupling to non- relativistic fermions [la], this strongly suggests 
that the putative deconfined phase in d = 3 is always unstable. Compact U(l) theories 
in d = 3, with or without matter, would appear therefore generically to be permanently 
confining. 

In the present article we study the issue of confinement in CQED3 using the variational 
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treatment of the anomalous sine-Gordon (ASG) theory, which is dual to the original cQED 3 . 
By working to the second order in fugacity and including the screening effects we find that 
monopoles are free at any effective temperature in the ASG theory (i. e. for any number of 
fermion flavours in cQED 3 ). This suggests that fermions are permanently in the confined 
phase, and provides an additional support to the renormalisation group results of Refs. 
and jig. 

We introduce the CQED3 theory and its dual sine-Gordon version in Section [H] In 
Section ITTT1 we discuss the lowest order variational calculation that neglects screening and 
point to its limitations. We then propose a generalised self-consistent approach that includes 
higher orders in monopole fugacity and allows for the screening effects in Section I1VI In 
Section |3 we present the calculations to the second order. A summary of our results is given 
in Section fVTl 



II. cQED 3 AND THE ANOMALOUS SINE-GORDON THEORY 

We will be interested in the phases of cQED 3 , with the gauge field coupled to massless 
relativistic fermions on a lattice: 

S[x, a] = S F [x, a] ~ 5-3 Yl cos ( F ^( x )) ■ W 

The sites of the three dimensional quadratic lattice are labeled by x M = {xi,x 2 ,r}. Here, 
is the usual field-strength tensor F^ u = A M a v — A^a^; the lattice derivative is defined 
by A At a J/ (x) = a l ,(pc + fl) — a u (x). Sp is the lattice action of massless fermions coupled to the 
gauge field which reduces in the continuum limit to QED 3 with Af flavours of four- component 
Dirac spinors. Using staggered fermions, this takes the form 

N { /2 

V 
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Sf[x, a] = jEE V, (x) Xn(x)e"v« X n(x + p) - Xn(x + /i)e-**«X„(x) (2) 



x,/i n =l 

where rji = 1, r]2 = (— l) Xl and 773 = (— \) X1+X2 jl7 |- 

In the case of continuum QED 3 , the fermion polarisation to one- loop order is 

n^(p) = §p{^-^}- (3) 

Incorporating compactness of a M in the spirit of Villain approximation Q], this suggests 
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that we consider a theory closely related to CQED3 




16 |A 



N f 



(F^(x) - 2im fJlU (x)) 



(4) 



where the n^ v are integers. The action (4) has the same continuum limit as CQED3 to the 
leading order in large iVf and may be understood as a compact quadratic approximation to 
it. In the remainder of the article we assume that the original cQED 3 and the theory (4) 
are in the same universality class. 

In the presence of fermions, when iVf 7^ 0, the original Maxwell term proportional to 1/cq 
becomes irrelevant at large distances, and can be neglected with respect to the second term 
in Eqn. (J£j). This action can be then be put into the alternative form (see Appendix IXJ) 



This is the partition function for a gas of monopoles of charge q a = ±1, interacting with 
a potential V^x). In our case, the potential has the form V(k) = 1/|^| 3 in Fourier space, 
which is the logarithmic interaction in three dimensions. 

The problem now appears to be rather similar to the two dimensional Coulomb gas, where 
the logarithmic interaction may result in the BKT vortex-antivortex binding transition. The 
mechanism of such a transition stems from a simple energy-entropy competition, as both 
entropy and the interaction energy are proportional to the logarithm of linear dimension of 
the system: at low enough temperatures, it is energetically favourable for opposite vortex 
charges to form bound pairs, while as temperature is increased, entropy at some point takes 
over, and vortex-antivortex pairs unbind. Although the form of the partition function (jJJJ) 
resembles that of the two dimensional Coulomb gas, it is by no means guaranteed that such 
a scenario will still hold in three dimensions. In particular, the effect of screening of other 
dipoles on the potential felt by two widely-separated monopoles, which is neglected in this 
naive energy-entropy argument, can drastically affect the result. 

To systematically address this issue, we first note that Eq. (JSJ) is equivalent to the partition 
function with the anomalous sine-Gordon (ASG) action (see Appendix lA"|) 



where the fictitious temperature is T = 2/(-7r 2 A r f), and y is the fugacity of the monopoles. 
The non-analytic gradient term proportional to |g| 3 is a consequence of the coupling of 
relativistic massless fermions to the gauge fields. 
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It is possible to construct an upper bound for Fa so, the free energy associated with the 
action using the Gibbs-Bogoliubov-Feynman |19( (GBF) inequality, which is discussed 
in the next section. We will argue that this self-consistent mean-field approximation to 
the free energy of the system unfortunately misses the screening effects of the medium, 
and consequently incorrectly suggests the BKT transition. An improved calculation that 
incorporates such effects is then formulated in the following section. 



III. VARIATIONAL APPROACH 

The GBF inequality imposes a strict upper bound on the free energy Fasg through the 
relation 

FASG < -Fvar = Fq + (SasG ~ S ) , (7) 

where Sasg is defined in Eqn. © and S is a trial action chosen to approximate Sasg; Fq is 
the free energy associated with So and (...) represents averaging within this ensemble. The 
trial action may be chosen to have the Gaussian form 

^o[0] = ^E^(q) G o 1 (q)0(-q), (8) 

q 



so it becomes particularly simple to calculate F var : 
F mr 1 f d 3 q 



(9) 



V 2 J (2tt) 3 

Minimising F var with respect to Go(q) yields the optimal Gaussian theory that approximates 

-^ASG : 

^ = 0^G 1 (q) = T|g| 3 + a, (10) 

with the 'mass' a determined self-consistently through 

1 f d 3 q 1 1 



a = 2y exp --/ 



2 J (2tt) 3 T|g| 3 + a 



Ik 

T 



A is the ultraviolet cutoff and T c = l/(127r 2 ). Determining the solutions of Eqn. (fTTj) 
amounts to identifying the roots of the function 

f(a) = a-2yll + — \ . (12) 



It is evident that a = is one such root for all values of T. We next demonstrate that a 
solution with finite a exists for T > T c . In the limit of small a, f(a) has the form 

fa, T <T C 

/(^« A H tit ° 13 
[ -a T ^ T , T>T C 

while for large a 

/(a»A 3 ) = a, VT. (14) 
For T > T c , /(a) changes sign and thus has a root with a > 0, while only the a = solution 



exists for T < T r 



2l|. 



The stability of the a = solution for T > T c can be determined from the variational 
free energy (jUJ) with the solution (fTUJ) for Gq . Evaluating the free energy we get 

= T C A 3 In (a + TA 3 ) -2y(l + ^ . (15) 

Then 

\ (F V M ~ F v U0)) = ° {Tc ~ T) + 0(a% (16) 

so that for T > T c any solution with a > is of lower free energy than with a = 0. That is, 
the stable solution at T > T c has finite a. 

To understand the physical meaning of the non-trivial solution it is useful to calculate 
the monopole density from the variational free energy 0: 

1 dF mr 



Pm 



V dp 

y dF va 



V dy 

(17) 

where we have used the definition of fugacity y = exp{/i}. We see that a is exactly the 
monopole density pm, so that o ^ may be identified with the plasma phase of free 
monopoles, while a = indicates the dipole phase. The simple variational calculation 
would therefore suggest that monopoles undergo a binding-unbinding transition at T = T c 
(i. e. at = A^ c = 24) in exact analogy with the equivalent calculation one can perform for 
the standard BKT transition. The value of T c also agrees with the simple energy-entropy 
argument that can be constructed for an isolated vortex jl^ . 
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An obvious objection to this simple calculation is that minimisation of the variational 
free energy (9) by construction cannot yield any momentum dependence of the self-energy, 
but can only determine its constant part, the 'mass' a. The renormalisation group [15| 
treatment of the ASG theory suffers from the same problem to the lowest order in fugacity, 
and would likewise naively suggest the BKT transition. The same holds for the direct 
perturbative evaluation of the self-energy in the ASG. However, it is easy to check that 
the self-energy does become momentum dependent to the second order in fugacity, with the 
leading analytic term ~ q 2 at low momenta. This is just what one would expect based 
on the simple electrostatic analysis of the problem where this term translates into the 
Coulombic interaction in real space (when y = 0). The presence of such a term would, 
however, drastically alter our present considerations. Indeed, if we add by hand the term 
Qq 2 with Q ^ in the denominator of the integrand in the self-consistent equation (jllj) . we 



for all T. Hence, the non-trivial solution would exist for all temperatures, exactly as in 
Polyakov's original treatment of the pure gauge theory. This is natural since Q ^ means 
that the original logarithmic interaction between monopoles is, even without free monopoles 
and only with a finite density of dipoles, screened into the Coulomb interaction for which 
the standard argument for the confined phase readily applies. 

In the next section we propose a modified self-consistent calculation which provides a 
systematic perturbative approximation to the free energy and which reduces to the GBF 
method to the lowest order. As we will see in Section|3 such an approach has the advantage 
of including the screening effects in a self-consistent way, therefore overcoming the limitations 
of the purely variational theory discussed in this section. 

IV. SELF-CONSISTENT PERTURBATIVE APPROACH 

There are many ways in which one may generalise the variational method of the previous 
section. For instance, one may add a second-order term — |((Sasg — Sq) 2 ) +\{S asg~ So)o to 
F var and extremise the new energy functional 19J. Such a second-order extension, however, 
has little variational justification. For a more systematic generalisation, we go back to the 



find 




(18) 
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GBF inequality (J7J) and exchange Sasg with So to find 

F< = F + (Sasg - So) < ^asg- (19) 
Extremising F< with respect to a quadratic action So yields 

(0(-q)0(q)) o =(0(-q)0(q)), (20) 

which is nothing but the equation for the exact propagator in the ASG theory. The right 
hand side (RHS) of the equation, on the other hand may be rewritten as 

<*HU«q)> - W "t-if AS>0 - < 21 > 

\ e /o 

with AS = Sasg — So. Eqn. (18) in this form may be understood as a self-consistent 
equation for the action S , which we may attempt to solve by expanding the RHS in powers 
of AS, for example. To the first order in AS this becomes 

(0(-q)0(q)AS) o - (0(-q)0(q)) o (AS) o = 0, (22) 

which is precisely the relation one would obtain from extremising F var with respect to So. 
That is, the first order approximation to Eqn. (J2U)) reproduces the GBF result from the 
previous section. 

Eqn. (|2*Uj) forms the basis of our modified variational approximation to -Fasg- To the first 
order in AS it reduces to the GBF equation of the previous section, and when solved self- 
consistently to all orders gives the best variational lower bound to the free energy, provided 
by F < in (fTH|) . In addition, consider the expansion of Eqn. (|2T)|) to order (AS) n . One can 
show (see Appendix IB 2|) that the resulting expression is the same as the one that would 
arise from extremising the function 

= 1 • (23) 

Here F^ stands for the expansion of the true free energy of the system, -Fasg ; i n powers 
of AS, truncated at (AS) n . Similarly denoting by F^ the truncated expansion of F< in 
Eqn. (fTUj) . it is not difficult (see Appendix IB 1|) to show 

T?{n) _ p( n ) 

Ffi = F (n) + —. (24) 

71/ 
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It is then clear that the sequence {Fj^} converges to Fasg f° r an Y So- Therefore, the S 
determined self-consistently from Eqn. ()2()j) also yields the variational sequence that best 
approximates Fasg within the family {F^fSo]}. 
To the second order Eqn. ()2U|) reads 

(0(-q)0(q)AS) c o - ^(0(-q)0(q)(AS) 2 ) C o = 0, (25) 

where both terms are connected averages given by: 

(0(-q)0(q)AS) c o = (0(-q)0(q)AS) o - (0(-q)0(q)) o (AS) o , (26) 
(0(-q)0(q)(A5) 2 ) C o (0(-q)0(q)(A5) 2 ) o -(0(-q)0(q)) o ((A5) 2 ) o 

-2(0(-q)0(q)A5) o (A5) o + 2(0(-q)0(q)) Q (A5) 2 . (27) 

We discuss the results of the second-order self-consistent Eqn. f)25|) for the ASG model (JBJ) 
in the next section. In particular, we will show that the density of free monopoles is finite 
at all T > 0, and that charge should consequently be permanently confined in cQED 3 . 



V. CONFINING SOLUTION FOR T > 



From the definitions of 5*asg an d So (Eqns. E]and|Bl) it is straightforward to calculate the 
connected averages of Eqns. (|2T)j) and (|2T|) . Our second order equation (|23j) then yields the 
quadratic equation for G$ x (q) 



where 



Go^q)! - A[q, G ]G \ q ) + B[k, G ] = 0, 



A[q,G ] = -T\q\ 3 + 3a + ab - 2a 2 c + ^(-l) n d n g 2n 

1 



B[q,G ] 

In Eqns. fl2T?l IHUj) . we have defined 



ZT 2 q 6 + 2aT\q\ 3 . 



a = ye 
b 

c 



4a>(o) 

H 3 k tT 1 



(I, 



J d 3 R[l -coshA)(R)], 

sinh D (R) 



■ (RcosP^ n 



d d R 



(2n)! 



(28) 

(29) 
(30) 

(31) 
(32) 

(33) 

(34) 



and the real-space propagator is D (H) = j d 3 k/(27r) 3 Go(k)e jk ' R . 

We can solve the quadratic of Eqn. (J2~%|) and expand in powers of |g| 3 /v4 to yield the 
result for Gq X (q): 

G x (q) = m + Q{m)q 2 + f(m)\q\ 3 + ■■■ (35) 
where the coefficients are defined as 

m = ^{A ±\A \}, (36) 
Q(m) = a 2 d, (l ± , (37) 

and with Aq = A[q = 0, Go]. For these equations, we should choose the solution correspond- 
ing to the upper sign in Eqns. (j3*o1 - I3*H )) to ensure that m > 0. In what follows, we neglect 
terms higher order in q than g 3 as they should be irrelevant at low momenta. 

As announced, the second order result includes additional renormalisation of the bare 
terms as well as the generation of new momentum dependent terms. Most importantly, the 
leading term proportional to q 2 has now appeared. 

In the analysis in Section ITTT1 we found that the bound phase of monopoles corresponded to 
low T. In what follows we will restrict ourselves to low temperatures by assuming TA -C Q 
and show that monopoles are unbound even for arbitrarily small temperatures. By continuity 
this would imply that they are free at all temperatures. 

Let us start by examining a: 

a = y exp j--£> (0) 



1 r d 3 k 1 



y exp 



2 J (2tt) 3 (Q(m)k 2 + m) 



=H-4^t A -^-n A v - ■))/■ (39> 

When m — > 0, we will assume m/Q(m) — > 0, and justify this assumption a posteriori. The 
coefficient a now takes the form 

"^"{"S^r}' ">«a 3 + o ( t). (40) 
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Next, we examine the equation for b 

d 3 k (T 



(*r)« V 2l*l 3 -^ -'(k))[G„(k)f 



-jA>(0) + 0(T). 



From this we find 



A 



m< A 3 + 0(T). 



47r 2 Q(m) 

Next, as the terms c and d§ always appear together, we consider the combination 



(c + 4) « | d 3 R M -expi-^ 



d 3 k 



,ik-R 



(2tt) 3 (Q(m)fc 2 + m) 



+ 0(T) 



Q(m) 



'1/' , )( , 



+ 0(T). 



Evaluating this yields 



(c + d ) = m-\ m<A 3 + 0(T). 
Similar analysis applies to the coefficient d\. 



1 



R 3 e~V 



dl = a / djR 

6 Jo Q{m) 



Q(m) 



+ 0(T), 



which gives 



Q(m) 



m < A 3 + 0(T). 



Evaluating the Eqn. (J37|) for Q then we find 



Q = 2a 2 dx 



2y exp 



A 1 Q 



2tt 2 Q m? 



+ 0{T). 



Solving this for Q ^ yields 



and we see that m/Q(m) indeed approaches zero as m — > 0, thus justifying our 
assumption. Substituting this solution for Q(m) into our mass equation gives 



m = A 
m 

* 71 



3 - v 7 ^ - In 



??? 
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which can finally be solved for m ^ to give the finite mass solution 

m * = V2e 2 ^~ 3 y. (50) 

The corresponding finite value of Q is 

Q* = (51) 

2tt 2 (3-2 v / 2) 

Note that m* is proportional to y so that small fugacity translates to small m*, in accord 
with our assumption that m-C A 3 . 

To show that monopoles are free when m / 0, we calculate the monopole density as in 
Section 11111 From Eqn. ()23|) we see that the free energy associated with our second order 
equation (j23|) is 

F% = F o+ (A5> - \(^S)\ + i(A5>;. (52) 
From this, the monopole density can be calculated. 

= 2a + ab — 2a 2 c. 



For m = 0, the monopole density vanishes, while for our finite m solution 

m* ( ,- 1 / 2m* \ 

2V2-1 



(2) 

Pm 



v^V 167rV(Q*)V 
> 0. (54) 



From the free energy (|52|1. it is also possible to show that the finite m solution is the stable 
solution for all temperatures. In fact, the free energy diverges as log(l/m) as m approaches 
zero, but has a finite value for finite m. It is then the free phase of monopoles which is 
favoured at all temperatures. 

Thus we have demonstrated that for arbitrarily low T a finite mass solution always exists 
for the self-consistent equations (|36l-l3*Sj). This implies that monopoles are always free at 
low temperatures, or, in terms of the original lattice model ([T|), that the electric charge is 
presumably confined for any number of fermion flavours. 



VI. CONCLUSION 



We have studied cQED 3 where massless relativistic fermions coupled to the compact gauge 
field result in a logarithmic interaction between magnetic monopoles. One may suspect 
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that this could lead to a BKT-like transition where free monopoles bind into monopole- 
antimonopole pairs at low enough effective temperatures. Although the simplest mean-field 
approximation would predict such a transition, we argued that by design this treatment 
misses the screening effects, argued to be crucial in this problem. To address this issue 
we developed a combined variational-perturbative approach which allowed us to include 
screening self-consistently. The modified theory then leads to the plasma phase of free 
monopoles as being stable at all temperatures, in agreement with the renormalisation group 
treatment of the problem 1 1 51 1 - 

n n 

CQED3 has been studied numerically in 228 an d .23]. Our calculation appears to be in 
agreement with the numerical results of (23], where only a single phase was observed. We 
hope that this and previous work on cQED 3 will motivate renewed efforts in this direction, 
using bigger system sizes that have recently become available 24 1 . 



VII. ACKNOWLEDGMENT 

This research has been supported by the NSERC of Canada and the Research Corpora- 
tion. 



APPENDIX A: MONOPOLE GAS AND THE SINE-GORDON ACTION 

In this appendix, we fill in the details in going from the action to the Coulomb gas 
partition function (jSJ) and the sine-Gordon action (JBJ). For this purpose, let us write (jU on 
the lattice in a more general form as 

S[a,n] = ^E(i^(x) - 2 7 rr V/ (x)) u (x,y) (F^(y) - 27rn^(y)), (Al) 

1 x,y 

where w _1 (x, y) = (16/iVf)|A x |5 X) y and A XiM / = /(x + ft) — /(x) denotes the lattice deriva- 
tive. Introducing an antisymmetric Hubbard- Stratonovich field M^ v we find 

S \ E M f( x ) y) M M + * E M^(x)(F^(x) - 2tt7v(x)) 

x,y x 

= jE^( x K 1 ( x ! y)^(y) + ^E(Vxa-27rn)(x)-b(x), (A2) 

x,y x 

with = t^xMyx and n M = e^n^x- Integrating over the gauge field constrains the b-field 
to be curl-free, so we can take it to be a gradient on the lattice b = A(p. Performing the 
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lattice version of integration by parts and integrating over <p yields 

S 7 E A /M x ) y) A MV(y) + 2« £ A ■ n ( x ) V?(x) (A3) 

x,y x 

-> ^E/'W^.y) p(y)' ( A4 ) 

x,y 

where v _1 (x, y) = — (l/87r 2 )A X)At tt _1 (x, y)A y ^ is the (inverse of the) potential and p = A M n M 
is the density of magnetic monopoles. Using the expression for w(x, y) we find, in the 
continuum limit, 

v^iVf r d 3 k e ik '( x - y ) 

v(x,y) = 



2 7 (2tt) 3 fc 3 



^(x-y). (A5) 



2 

Thus, for a system of iV monopoles with a density p(x) = Y,a=i QaS(x — x a ), we obtain 

Smon = X <laqbV(Xa ~ Xft) , (A6) 

a, b 

as in Eqn. (jSJ). 

We now proceed to show that (jA6|) with q a = ±1 is equivalent to the sine-Gordon 
action (j§J. Higher charges are irrelevant for large enough Nf. To this end, let us introduce 
the bare action 

S h = i | d 3 x d 3 y ^(x - y)0(y), (A7) 

so that u(x — y) = (0(x)0(y))b- Showing the fugacity by y, we may write the monopole 
partition function in the grand-canonical ensemble as 



7 

^mon 



y N r N f l 1 

J2j^ JU d3x i E ex p\ - 2E»W x »)^W>b \ 

N • i {q a =±l} { a,b ) 



E^f/lK* E exp 

N JV - J i {q a =±l} \ 

exp <.2y J d 3 x cos0(x 



E<?a 



x„ 



= Z^^asg, (A8) 

where Z h is independent of (j> and 

Zasg = y X>0 exp {-- J d 3 x d 3 y ^tT^x - y)0(y) + 2y J d 3 xcos0(x)| . (A9) 

Inserting the definition of ?; _1 (x — y) into this last expression we immediately arrive at the 
anomalous sine-Gordon action ©. 
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APPENDIX B 



1. F< n) 



Here we will show that indeed satisfies Eqn. (124)) . To this end, first let us define 
AF^' = F( n+1 ) — F^ n \ Then we may equivalently show that 

= F + AF (1) + 2AF^ + ■■■ + nAF (n) . (Bl) 

Let us also denote the path integral over the field 0(q) by Tr and define, for a real variable t, 

T{t) = - In Tr (e~ s °e- tAS ) . (B2) 



Then, J-"(l) = — InTr exp(— S) = -Fasg an d 

dJF(t) Tr (AS e~ s 



dt 



(AS). (B3) 



Tr (e-S) 

On the other hand, we may expand the RHS of Eqn. (|B2ft in powers of AS as 



r(t) = F + y £AJ* i \t) (B4) 

i=i 



where A^(t) = f A^W(l) = f AF«. Thus 

&T(t) 



(AS) 



dt 



f^iAF^. (B5) 



i=l i=l 

Upon insertion of Eqn. ()B5|) into the definition of F< in Eqn. f)19|) and truncating the 
expansion at % = n we find (JET]) . 

2 F (n) 

^ • -T var 

In this appendix, we will give the proof for our claim that the extremum of F^} as defined 
in ()23|) is given by the expansion of Eqn. (|2(J|) to order (AS) n , i. e. 

= <=> (0(-k)0(k)) o = (0(-k)0(k))W. (B6) 

The calculations are, for general n, cumbersome and not very instructive so we will first 
present the case for n = 2 which is also the one with which we are concerned in Section IVl 
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Setting n = 2, we see that Eqn. (J23j) is readily found by an expansion of the RHS of 
Eqn. (}2"T]) . To show that the same result arises from extremising F§v, it is first useful to 
establish 

= (sm) = 2w» (0(_q)0(q))o ' (B7) 

5 (9)o _ 5F . ,/ Sg 6S x 



<5G (q) 5Go(q)^ )o+ \5G (q) ^G (q)/ ' ^ 
where g = g(So) is an arbitrary function of Sq. Thus choosing appropriate forms of g for 
F« =F + (AS) and F( 2 ) = _ ±((AS) 2 > + |(AS)J we find 



5G (q) 5G (q) 

^(^l) (AS)2 > -(^) AS > (AS> " (mo) 

Inserting Eqn. ()B7|) into Eqns. (jB9l IB10|) and adding them we find that the restriction 
5F^/5Go(q) = leads to the same equation as Eqn. (|2*5Jl . 

The proof for arbitrary n goes along essentially the same steps as above. Various trun- 
cated expansions we have defined can be read off the Taylor expansion identity 

_ hMe ^^ + ^^ r (M, (B11) 

1=1 1 = 1 {k a } 



by setting i to the desired order. In (jBTTj) St, and V give an arbitrary splitting of the action 
into a bare and potential part respectively and 

{fc Q } fei=l fc;=l 

Notice that in Eqns. ()B9|) and (|B10|) all the terms are to the same order of AS, which is 
also the largest in the corresponding expansion of the free energy. By choosing Sb = Sq and 
V = AS and setting i = n in (IB 11 J) one can see, after some lengthy algebra, that the same 
is true for arbitrary n: 

^ {n) = A, ]n+l ^> ((AS) k %---((AS) k % 6F 
^o(q) ti { t} hl-'-hl 5G (cO 

_ A ir+ . r ^^fo(^) fcl )o(( A ^ 2 )o---(( A ^)o (m3) 
h { t} fci!^!---^! 
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Let us now define, for a real variable t, 

0(k,t) = _i n Tr{ e - s °- A5 -^ ( - k) * (k) }, (B14) 

so that dg(k,t)/dt\ t=0 = (0(-k)0(k)). Then, taking S h = S and 1/ = AS + t0(-k)0(k) 
in Eqn. (jBll|) to compute this derivative, it can be shown through additional tedious but 
straightforward algebra that 

(0(-q)0(q))W - (0(-q)0(q)) o = -2^2^ ( q )] 2 —^- (B15) 

where we have also made use of Eqn. (jB13|) . Thus, the requirement that F^) be an extremum 
implies Eqn. (|20p truncated at nth order, and vice versa, proving our claim (jB6|) . 
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